>> im=imread('etest.bmp'); %Bild einlesen, etest.bmp ist ein Graustufenbild
>> Imed=medfilt2(im); %Medianfilter
>> [rows, cols] = size(Imed);
>> Imed=double(Imed);
>> dx=[ Imed(1:rows, 1:cols) ]-[ Imed( 1:rows,1) Imed(:,1:cols-1) ]; % Erste partielle Ableitung nach x
>> dy=[ Imed(1:rows, 1:cols) ]-[ Imed( 1,1:cols) ; Imed(1:rows-1,:) ]; % Erste partielle Ableitung nach y
>> dxx=[ dx(1:rows, 1:cols)]-[ dx( 1:rows,1) dx(:,1:cols-1) ]; % Zweite partielle Ableitung nach x
>> dxy=[ dx(1:rows, 1:cols)]-[ dx( 1,1:cols) ; dx(1:rows-1,:) ]; % Gemischt partielle Ableitung von x nach y
>> dyx=[ dy(1:rows, 1:cols)]-[ dy( 1:rows,1) dy(:,1:cols-1) ]; % Gemischt partielle Ableitung von y nach x
>> dyy=[ dy(1:rows, 1:cols)]-[ dy( 1,1:cols) ; dy(1:rows-1,:) ]; % Zweite partielle Ableitung nach y
>> A=cell(rows,cols);
>> B=cell(rows,cols);
>> for i=1:rows,
for j=1:cols,
H=[dxx(i,j) dxy(i,j);dyx(i,j),dyy(i,j)]; % Hesse-Matrix
W=eig(H); % Die Eigenwerte der Hesse-Matrix
A{i,j}=W;
[V,D]=eig(H);
B{i,j}=V; % Die Eigenvektoren der Hesse-Matrix
end
end